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Abstract 

The procedure to improve the convergence in transverse momentum 
space of the NLL BFKL kernel using a tj-shift is revisited. An accu- 
rate approximation to this shift only depending on transverse momenta is 
presented. This approximation is based on a Bessel function of the first 
kind with argument depending on the strong coupling and a double log- 
arithm of the ratio of transverse scales. A comparison between different 
renormalization schemes is also included. 

1 Introduction 

In recent years there has been an intense activity trying to understand which are 
the effective degrees of freedom driving the strong interaction at large energies. 
For processes where the center-of-mass energy is much larger than the other 
scales in the scattering an instructive picture emerges from the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) approach pQ. The effective picture in this framework 
relies upon f-channel "Reggeized" gluons interacting with each other via stan- 
dard gluons in the s-channel. Although this simple structure must be modified 
at higher energies in order to introduce unitarization corrections, there is a win- 
dow at present and future colliders where the BFKL predictions should hold. 

The leading logarithmic (LL) approximation in BFKL rcsums terms of the 
form (a s lns) ra to all orders. In such a limit diagrams contributing to the run- 
ning of the strong coupling do not appear and a s remains as a constant free 
parameter. Also the factor needed to scale the energy in the logarithms is not 
fixed, therefore the predictability of the LL approximation is limited. The situa- 
tion improves already in the next-to-lcading logarithmic (NLL) approximation 

'Alexander von Humboldt Research Fellow 



where diagrams with an extra power in the coupling without introducing an 
extra logarithm in energy are taken into account. In this case the coupling is 
allowed to run and the energy scale is fixed. 

The calculation of the NLL corrections to the BFKL equation was lengthy 
and completed in Ref. It was needed to evaluate corrections to the gluon 
trajectory [3] and real emission kernel 0], as well as to check that the boot- 
strap conditions are fulfilled in order to ensure that the complete set of NLL 
contributions are indeed correctly resummed 

Lately, to be able to perform phenomenological studies, the NLL calculation 
and numerical implementation of particular impact factors, describing the na- 
ture of the interaction between the Reggeized gluons and the scattering particles 
for different physical processes, is a remaining challenge [H]. 

In the construction of BFKL cross-sections the gluon Green's function is 
the element carrying the energy dependence and therefore it deserves a careful 
study. There has been a lot of theoretical work trying to understand its structure 
which, since the appearance of the NLL kernel, has revealed as quite non- 
trivial. It was shown in Ref .7 that, in the transverse momentum space, the 
NLL Green's function presents some instability. The origin of this instability 
and more convergent kernels including renormalization-group (RG) generated 
terms were presented in Ref. Based on this type of RG-improved kernels, 
including more sophisticated approaches, several studies of the Green's function 
have been performed in Ref. . Other analysis of different aspects of the NLL 
kernel can be found in Ref. |10| . 

The aim of the present work is to revisit the approach of Ref. [H] and to 
extract the structure in transverse momentum space of the double logarithms 
there resummed. In order to do so, an approximation to the original w-shift 
will be performed. It will be shown how this approximation is an accurate one 
and the corresponding expression in transverse momentum space very simple. 
The main feature of the new RG-improved kernel will be that it does not mix 
transverse with longitudinal momentum components. From the practical point 
of view it also allows for its straightforward implementation in the method 
of solution of the NLL BFKL equation proposed in Ref. (for reviews see 
Ref. [12). 

In Section|21the NLL BFKL kernel will be presented and its representation in 
7-space will be introduced. The particular RG-based scheme under investiga- 
tion will be described at the end of this section. In Section OH the approximated 
solution to this scheme will be derived and an analysis of its accuracy in the 
minimal substraction (MS) renormalization scheme will be shown. In Section^] 
the corresponding RG-improved kernel is derived in transverse momentum space 
and its asymptotic behaviour in terms of the ratio of transverse scales discussed. 
Section [S] is a brief one analysing the effect of changing the renormalization to 
the gluon-bremsstrahlung (GB) scheme. Finally, the Conclusions are presented. 
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2 NLL BFKL kernel and modifications based on 
the renormalization group 

The starting point of this analysis is the following representation of the NLL 
BFKL kernel acting on a smooth function 



J d 2 ^ 2 IC(qi,q 2 ) f (ql) 
d 2 q 2 



(1) 



s 



f («§) - 2- 



ft 2 




Q1-Q2 



^T(^f) + ln*(| 



fill) 



where a s = g 2 ^/^ 2 ), ft = (11N C - 2n f ) /3, 5 = (4 - tt 2 + 5/3 /iV c ) /12, 
and the function T(q 2 , q 2 ) can be obtained from Ref. |2j- As it stands this ex- 
pression is written in the MS renormalisation scheme. Alternatively, as in other 
calculations where a resummation of soft gluons is involved, the GB scheme ^3] 
can also be used: 



a 



GB 



a s (1 +Sa s ). 



(2) 



Both schemes will be under analysis in this work. 

As it is well-known the collinear structure of this kernel can be extracted by 
projecting it on the modified LL eigenfunctions (a s (g 2 ) /a s ([i 2 )) 
i.e. 



d 2 q 2 lC(qi,q2) 



(ll) 



as (ql) 



4) 7 =«.(fl?)xo(7) + ^Xi(7). (3) 



Note that the prefactor (a 8 (q 2 ) l a s(q 2 )) 



-1/2 



(q 2 ) (f3 /8N c )\n(q 2 /q 2 ] 



generates a NLL term which makes the 7-representation of the kernel invariant 
under the 7^1 — 7 transformation . 

Although a s (<? 2 ) ~ a s (^ 2 ) ^1 — a s (/i 2 ) j^- In breaks the scale invari- 
ance, for the collinear analysis of the kernel it is enough to work with the scale 
invariant pieces which, at LL and NLL accuracy, read 



Xo (7) 



2V(l)-V(7)-^(l-7), 



(4) 



Xi (7) = S\o (7) + \ O" (7) + (1 - 7)) - ^ 



•(7) + 0(1 -7)) 



7r 2 cos (Try) 
4sin 2 (7T7)(l - 27) 



(5) 



1 + ^ (2 + 37(1 - 7) M + -Cs - ^Xl (7) 
JV C 3 J (3-2 7 )(l + 2 7 )j 2^ 3 8N c Xo(1) 
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All these expressions are re-written here just to set the notation, ip (7) = 
V (7) /r (7) and 

0(7) + 0(1 -7) = 



m=0 



7 + to 1^7" 



For future reference it is important to remark that the terms with second deriva- 
tives in the ip function directly stem from the double logarithm in Eq. . The 
pole structure of this kernel around 7 — 0, 1 is as follows: 

Xo(7) - - + {7^1-7}, (7) 

7 

cL b 1 , 

Xi (7) * - + ^-2^ + {7-l-7}. (8) 

The cubic poles come directly from ip" and 

5 /3q 13 n f 55 1 0a rif 11 



12iV r 36 N? 36' 8N r 6N? 12' 



(9) 



Throughout this paper for numerical analysis the value n/ = 3 will be taken. 

In this 7-space representation it is known that the cubic poles do compen- 
sate for the equivalent terms appearing when the symmetric Regge-like energy 
scale so = <?i<72 is shifted to the so = q\ 2 choice in DIS-like processes. This 
compensation takes place at NLL, the accuracy at which the BFKL kernel is 
known. It turns out that higher order terms beyond NLL, not compatible with 
RG evolution, are also generated by this change of scale. The NLL truncation 
of the perturbative expansion is then the reason why the gluon Green's function 
in the BFKL formalism develops oscillations in the q\jq\ ratio when this ratio 
is very far from unity. Along these oscillations the Green's function can have 
negative values. 

At present collider energies it remains to be seen how important these os- 
cillations are when a full NLL BFKL-resummed calculation of a physical cross 
section, including impact factors, is carried out. Nevertheless there have been 
attempts in the literature to improve the behaviour of the BFKL resummation 
in the q\jq\ variable. One of the original proposals, the subject of the present 
analysis, shows how it is possible to remove the most dominant poles in 7-space 
incompatible with RG evolution by simply shifting the w-pole present in the 
BFKL scale invariant eigenfunction. This shift has to be performed with care 
not to double count terms and to match the NLL accuracy of the original BFKL 
calculation. In the present study the focus will be on the scheme proposed in 
Ref. 0: 



LO = 

a s (l + U + y ) a^j (2^(1) - V- (7 + I - ba s ^ 



.) -V(l-7+|-ba s )) 



( \.| (7) ! ( i Xo ( 7 )-b) (^,'(7) +<//(!- 7))- (a +y ) \„h>) . 
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By iterating this expression new terms beyond the original BFKL calculation 
appear improving the convergence in q\lq\ of the expansion. It is worth noting 
that the shift immediately resums these new terms to all orders in the coupling. 
Therefore, any attempt to modify this approach should retain this feature. 



3 "All— poles" resummation 

The main idea underlying the present analysis is that the numerical solution to 
the equation 

lo = a s (2^(1) -^(7+ |) ~i> (l-7 + I)) (11) 
can be approximated remarkably well by the all-orders resummed expression 

w = /t#s{ ( x7_1 + x ~ 7 ) {Wx Jl ( ^ ln2 x ) ~ 2as } ' (12) 

where Ji(z) is the Bessel function of the first kind. 

Before showing the calculations leading to this formula the expression in 
Eq. H12|) is compared to the numerical solution of Eq. in Fig.^ The values 
corresponding to the LL BFKL kernel are also included in the plot. The w-shift 
removes the poles at 7 = 0, 1 and affects the saddle point region by reducing the 
LL intercept from about 0.55 at a fixed coupling of a s = 0.2 to a value close to 
0.31. The approximated value obtained from the expression containing Ji(z) is 
0.35. It will be shown below how this approximation to the shift is even more 
accurate when the full NLL corrections are included. On the last plot of Fig. ^ 
one can see that the approximation is a stable one under variations of the strong 
coupling values. 

In the following the full NLL scale invariant terms are taken into account. 
In general one can consider the w-shift of the form 

uj = a s {l + Aa s ) (2^(l)-^( 7 +! + Ba s ) - i/j (l - 7 + | + B a.)) . (13) 

It is convenient to use the following representation of the LL eigenvalue of the 
BFKL kernel in Eq. Q: 

»(7) = Ef^ + r^r xt)- ( 14 ) 

* — ' \7 + m 1 — 7 + m m+1 J 

rn—O v ' ' ' 

This representation reveals the singularity structure of Xo(7) which has an in- 
finite number of poles at the 7 = —to, 1 + to points along the real axis. In this 
way the general shift can be written as 



a s (1 + Aa s ) 



(15) 

1 1 2 \ 



7 + m+|+Ba s l-7 + m+f +Ba s m+1. 
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Figure 1 : Behaviour of the 7-representation of the LL BFKL kernel compared 
to the w-shift of Eq. 111(1 . For a fixed coupling of a s = 0.2 on the top and 
middle plots around the saddle point, and at the dominant region 7 = 1/2 for 
different perturbative values of the coupling (bottom). The approximation to 
the shift using the Bessel function resummation as in Eq. is also shown in 
the plots. 
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In those regions close to the 7 ~ —m poles a good approximation to Eq. (|15fl is 

a s (1 + Aa s ) 



7 + m + I + B a s 

with the straightforward solution 



(16) 



u = -( 7 + m + Ba,) + l7 + m + Ba s l f 1+ 2& »( 1 + A& >) \ . (17) 

\ (7 + rn + Ba s )'y 

The same logic applies to the regions around the 7 ~ 1 + m poles. In the 
following it will be shown how, to a very good accuracy, the solution to the shift 
in Eq. Q15[) can be obtained by simply adding all the approximated solutions 
at the different poles plus a term related to the virtual contributions of the LL 
BFKL kernel, i.e. 



ui = 

m=0 



£ \ -(l + 2m + 2Ba s ) + | 7 + m + Ba s | ( 1 + 2 " s (1 + A '° 



(7 + m + B a s ) 2 



, 2a s (l + Aa s ) V 2a s (l + Aa s ) ■ 

|l- 7 + m + Ba.| 1 + - — — -—2 — ^.(18) 

(1 — j + m + Ba s ) I m + i 



From now on this expression will be indicated as the "all-poles" approximation. 
In the region of interest, < 7 < 1, Eq. (|18() simplifies because 

|7 + m + Ba s | = 7 + m + Ba s , (19) 
|l-7 + m + Ba s | = l-7 + m + Ba s , (20) 

and one can also introduce the all-orders expansion of the square roots: 

2a s (l + Aa s ) 



E 



m + 1 

, n+l 



^ (-l)"(2n)! (as+Aal) f 1 . 

+ g 2«„! (n + I)- (7 + m + Ba.r + {7 " 7} J J ■ (21) 

At the 7 = 0,1 poles this expansion generates the NLL terms: 

a s o / A B 1 \ , . , . 

w~-L+a?t +{7^1-7}- 22 

7 V7 7 2 7 J 7 

Therefore, to match the original kernel at NLL, it is needed to set A — a and 
B = -b from Eq. ©. 

To demonstrate that this expansion improves the convergence of the BFKL 
calculation in the same way as the original w-shift, it is needed to include the 
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full NLL scale invariant kernel without double counting terms, i.e. 



u 



(23) 



{OO / oc 

m=0 L \n=0 





\7 + m (7 + m) 2 2(7 + m) 3 



) 



+ {7^1-7} 



7 + m 



In this way, as in the original w-shift, the NLL results are untouched and only 
higher order terms are added to the kernel. The extra terms being resummed 
to all orders in a s will be important to find a closed form for the kernel with 
corresponding eigenvalues as in Eq. (|23|l . 

It turns out that when the matching to NLL is performed the "all-poles" 
result reproduces the w-shift very accurately. This can be seen in Fig. El where 
the full ui— shift and the "all-poles" kernel are compared, together with the LL 
and NLL results which are also included. The fact that the imaginary part of 7 
at the maximum of the NLL scale invariant eigenvalue (middle plot of Fig. [2]) is 
not zero results in the oscillations in the qf/q 2 variable. This feature is removed 
when the RG-improved kernel is used, as it also happens for the "all-poles" 
kernel. The estimated intercept from the w-shift is of about 0.27, while that 
coming from the "all-poles" kernel is 0.26 (for 6t s = 0.2). The similarity between 
the "all-poles" kernel and the original w— shift remains for different values of the 
coupling (see bottom plot of Fig.EJ. 

4 Bessel representation in transverse momen- 
tum space 

The most striking feature of Eq. (|23|) is that the w-space is decoupled from the 
7-representation, i.e. they no longer mix as it occurs in the original w-shift. 
This has the implication that it should be possible to find an expression for 
a RG-improved BFKL kernel which does not mix longitudinal with transverse 
degrees of freedom. To find such a kernel it is convenient to define the following 
function related to a simple representation of the LL kernel: 




(24) 



2^(1) - V(7 " a„b) - V(l - 7 - u s b) 
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Figure 2: The 7-representation of the LL and NLL scale invariant kernels show- 
ing the instable behaviour of the last one. The RG-improved kernel by a w-shift 
is also included, together with the "all-poles" approximation proposed in the 
text. 
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Keeping in mind the coefficients of the expansion in Eq. (|23[l the following 
summation can be constructed: 



E 



(-l)"(q s + aq 2 )" +1 d 2n 
2 n n!(n + l)! d"f 2n 



n( 7 ) 



(25) 



dq 2 



\k 2 -q 2 \ \k 2 



2 \ t- 1 - ha !> L k^T °° 



E 



(-l)»(5 s +aa 2 ) " +l 
2 n n\ (n+l)! 

. n+l 



= ff (-1)"(2»)! (°. + adj) 

^^2 I1 n!( Ji + l)!(7 + m-ba s ) 2 "+ 1 11 ,J ' 

m— n— 1 v / \ ' / 



In order to include the n = case it is needed to add a term corresponding 
to the virtual contributions in the LL BFKL kernel, ensuring, in this way, the 
infrared finiteness of the final result, i.e. 



dq 2 



\k 2 -q 2 \ 




-2 (a. + &ai 



i i_ - \h — -q\ 

2 \ 



min(fc 2 , q 2 
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k 2 
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m— 
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a a; 



771+1 
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— ' ^ 2™n! (n+l)! (7 + rn-ba s ) 2n + 



(26) 



By introducing this extra n — piece it is then possible to resum the expansion 
in terms of a Bessel function and Eq. i|2ti|) is equivalent to 



dq 2 



-2 (o! s + a a 



2 s min(fc 2 , g 2 ) 
fc 2 + g 2 



(27) 



2 X 7-l-bfi s ^|i 



\ 



2(a s 



1^ | 



U L/2(a s + aa 2 )ln 2 



In order to match this expression to the NLL BFKL formulation it is sufficient 
to perturbatively expand Eq. (|2"7|l up to order a 2 : 



2 (a s + a a 2 ) 



Ji[ J2(a s +aa 2 )ln 2 (|2 



(28) 



a s + aa. 



v. -2 l fc _ l\ , 

ba 2 ^ Ljn 



It is important to note that this expression does not depend on the angle between 
the q and k two-dimensional vectors, therefore this scale invariant resummation 
can be directly added to the original NLL BFKL kernel without having to 
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angular average. The "all-poles" approximation only affects the zero conformal 
spin sector of the theory. This is particularly interesting for the study of final 
states where angular correlations are relevant. Moreover, these new higher order 
terms do not affect the cancellation of infrared divergences in the original kernel. 

It is also noteworthy to point out how the perturbative expansion of the 
"all-poles" kernel naturally generates the double logarithm already present in 
the original BFKL kernel (see Eq. QJ). What this resummation then achieves 
(as in the original w-shift) is to resum to all orders those most dominant trans- 
verse logarithms which will make the expansion compatible with renormaliza- 
tion group evolution in the collinear and anticollincar limits relevant in DIS-like 
configurations. 

For the sake of clarity, the only modification needed in the full NLL kernel 
to introduce the "all-poles" resummation is to remove the term 

a? 1 W£) (»> 



4 (q-k) 2 \k 2 
in the real emission kernel, IC r (^q, kj , and replace it with 




-ha 



2(a s + aa 2 ) 2 f q 2 

Ji a 2 (a s +aa|)ln 1 



k 2 

an; fb^^lnf^L ,:-!,„ 



From the asymptotic behaviour of the Bessel function, when the coupling is 
small and the difference between the q 2 and k 2 scales is not very large then 



and the influence of the "all-poles" resummation is minimal. Therefore the 
"Regge-like" region in momentum space is not largely affected. But when the 
logarithm in the ratio of transverse momenta becomes larger due to a larger dif- 
ference between the scales entering the real production vertex, then the asymp- 
totic behaviour develops the oscillatory form 




Ji 



automatically compensating for the unphysical oscillations present in the origi- 
nal formulation of the NLL BFKL equation. 





11 



5 Change of renormalization scheme 



In this section the study of the effect of changing the renormalisation point, and 
using the gluon-bremsstrahlung (GB) instead of the MS scheme, is performed. 
This is equivalent to the following redefinition of the position of the Landau 
pole: 

A G b = A m exp(s^J. (33) 

With this choice the coefficient for the double poles in the NLL kernel does not 
change while the one for the simple poles now reads 

7T 2 13 71/ 67 
a " 12~36iVf _ 36- (34) 

In Fig. |2| this GB scheme is under analysis. There it is shown how again the 
"all-poles" kernel is a very good approximation to the w-shift. With a MS 
value of a s = 0.2 the corresponding coupling in GB scheme is ~ 0.23. The 
value of the intercept at LL is of ~ 0.64 while the RG-improved prediction is 
of ~ 0.29, which coincides with that steming from the "all-poles" analysis. It 
is clear that the influence of the choice of renormalisation scheme will be larger 
when the full kernel is introduced and running coupling effects are also taken 
into account. 

As a final remark, if running coupling effects are to be included, it is inter- 
esting to note that Eq. corresponds to an angular-averaged representation 
of the gluon Regge trajectory and the real emission kernel. If this angular in- 
tegration is not performed then the cancellation of infrared divergencies can be 
written as 



-' (( ? 2 ' A ) + / ^r cusp (a s (fc 2 ))0(fc 2 - A 2 ), (35) 



where the gluon Regge trajectory in this regularization reads 

ri 2 dp 

w (g 2 ,A) = -J ^ —r cusp (a s (fc 2 )) + constant. (36) 

The constant of integration can be fixed to be a 2 from the first reference 
in [2] and 

r cusp (a s (fc 2 )) = a s (k 2 ) + a 2 s S, (37) 

is the so-called "cusp anomalous dimension" of Ref. 14 . In principle it is 
possible to use a s (k 2 ) = 4iV c /(/3o In fc 2 /A 2 ), multiply the NLL equation by the 
extra logarithm and obtain the gluon Green's function as the solution to a 
diferential equation in 7-representation. An alternative is to solve the equation 
as proposed in Ref. |llj . 
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One of the targets of the present paper is to show how it is possible to intro- 
duce a RG-improved kernel which still allows for an exact solution of the NLL 
BFKL equation as that proposed in Ref. ^J. The solution there calculated is 
based on the iteration of the BFKL equation for the partial wave expansion of 
the Green's function. This iteration has a multiple pole structure in the complex 
w-space which allows for a simple inversion back into energy space. This inver- 
sion using a Mellin transform would not be possible if an extra w-dependence in 
the kernel, as the one in the w-shift, is present. In this context the prescription 
proposed here in Eq. H29|) and Eq. (|3L)[) allows for a straightforward implemen- 
tation of the RG improvements given that it just modifies the real emission 
part of the NLL BFKL kernel without introducing any extra w-dependence. 
Therefore the iterative procedure used in Ref. remains valid. In particular, 
if a different choice for the argument in the running coupling is considered the 
method used in Ref. ^5 always provides the exact solution. The study of all 
these effects at the level of the gluon Green's function will be performed in a 
coming publication. 

6 Conclusions 

The original analysis of Ref. [H] for the improvement of the convergence in 
transverse momentum space of the NLL BFKL kernel is revisited. An intrin- 
sic feature of such renormalization group improved kernels is that the variables 
Mellin-conjugated of the energy and the transverse momenta, uj and 7 respec- 
tively, do mix due to the so-called w-shift. This makes the study of these 
schemes quite complicated. In the present work it is shown how it is possible to 
avoid the difficulties of the w-shift without loosing the physical insight. This is 
done by approximating the solution to the original equation at each of the poles 
of the original kernel. The result in 7-representation is that of Eq. (|23() named 
as the "all-poles" resummation. It turns out that the corresponding scale in- 
variant kernel in transverse momentum space with eigenvalue as in Eq. I|23|) 
can be found and has a very simple form. Its structure in terms of a Bessel 
function of the first kind is shown in Eq. (|3U|) . The main feature of this kernel 
is that it only depends on transverse momenta and not on longitudinal degrees 
of freedom. 

This representation of the collinear improved kernels can be implemented 
immediately using the iterative approach described in Ref. m . The advantage 
of this approach is that it takes into account all physical effects, including the 
running of the coupling, exactly. It also allows for the study of final states and 
angular correlations. 
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